ESD-TR-72-30 


ti 

n» 


% 


Technical  Note 


1972-6 


R.  D.  Ruquist 


Electromagnetic  Scattering 
from  Turbulent  Boundary  Layers* 


tv 


K'lTUa'x  i  1 


JENTJFIC  &  TECHNICAL  ifcrOCMATlON  DIVISION 


21  January  1972 


Prepared  for  the  Advanced  Research  Projects  Agency 
under  Electronic  Systems  Division  Contract  F19628-70-C-0230  by 

Lincoln  Laboratory 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
Lexington,  Massachusetts 


f 


i 


AD  7  #612- 


Approved  for  public  release;  distribution  unlimited. 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 


LINCOLN  LABORATORY 


ELECTROMAGNETIC  SCATTERING 
FROM  TURBULENT  BOUNDARY  LAYERS 


R.  D.  RUQUIST 
Group  35 


TECHNICAL  NOTE  1972-6 

21  JANUARY  1972 


Approved  for  public  release;  distribution  unlimited. 


LEXINGTON 


MASSACHUSETTS 


The  work  reported  in  this  document  was  performed  at  Lincoln 
Laboratory,  a  center  for  research  operated  by  Massachusetts 
Institute  of  Technology.  This  work  was  sponsored  by  the 
Advanced  Research  Projects  Agency  of  the  Department  of 
Defense  under  Air  Force  Contract  F19628-70-C-0230  (ARPA 
Order  600). 


ABSTRACT 


The  incoherent  scattering  cross  section  of  an  underdense 
turbulent  boundary  layer  adjacent  to  both  a  perfectly  conducting 
ground  plane  and  a  perfectly  conducting  cylindrical  surface  is  de¬ 
termined.  It  is  found  that  the  effect  of  the  metallic  surfaces  on 
scattering  from  the  turbulent  plasma  depends  on  the  thickness  of 
the  boundary  layer  and  the  characteristics  of  the  incident  electro¬ 
magnetic  waves.  For  optically  thin  boundary  layers,  a  plane 
metallic  surface  will  enhance  scattering  in  one  linear  polariza¬ 
tion  by  up  to  12  dB  while  eliminating  scattering  in  the  other 
polarization. 


Accepted  for  the  Air  Force 

Joseph  R.  Waterman,  Lt.  Col. ,  USAF 

Chief,  Lincoln  Laboratory  Project  Office 


» 


I 


CONTENTS 


Abstract  iii 

I.  Introduction  1 

II.  Scattering  from  Turbulent  Media  in  the  Geometrical 

Optics  Limit  2 

III.  Scattering  from  Turbulent  Boundary  Layers 

on  Ground  Planes  5 

IV.  Scattering  from  Turbulent  Boundary  Layers 

on  Perfectly  Conducting  Cylinders  8 

Appendix  1 1 


v 


4 


I 


9 


I 


ELECTROMAGNETIC  SCATTERING 
FROM  TURBULENT  BOUNDARY  LAYERS 


I.  INTRODUCTION 

Scattering  of  electromagnetic  (EM)  waves  from  plasmas  residing  in 
turbulent  media  has  largely  been  treated  as  a  scalar  problem  rather  than  a 
vector  problem.  In  the  vector  problem,  the  radar  cross  section  (RCS)  of  a 
turbulent  medium  is  determined  from  a  coherent  summation  over  the  rays 
scattered  from  the  medium,  whereas  in  the  scalar  problem,  the  phase  of  the 
scattered  ray  is  ignored  and  an  incoherent  summation  is  used.  Generally, 
polarization  is  also  ignored  in  the  scalar  problem  because  the  incoherent 
summation  makes  all  polarizations  equivalent. 

Scattering  of  EM  waves  from  wakes  of  hypersonic  vehicles  has  been 
treated  as  a  scalar  problem  because  the  turbulent  fluctuations  in  the  "thick" 
wake  rapidly  randomize  the  phase  of  any  signal  propagating  through  the  me¬ 
dium.  However,  the  thickness  of  turbulent  boundary  layers  on  hypersonic 
vehicles  is  often  much  smaller  than  the  wavelength  of  the  incident  EM  waves 
and  one  would  not  expect  phase  to  be  randomized.  Therefore,  a  vector  solu¬ 
tion  is  appropriate. 
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Yakimenko  has  suggested  a  very  convenient  relationship  for  obtaining 
such  vector  solutions.  He  applied  his  technique  to  "thick"  turbulent  media 
and  discontinuous  media.  A  generalized  form  of  his  relationship  derived  in 
the  Appendix  may  be  used  to  investigate  "thick"  as  well  as  "thin"  media.  In 
Section  II  we  show  that  vector  solutions  are  essentially  limited  to  cases  where 
the  Born  Approximation  is  valid,  and  that  the  Born  Approximation  always  re- 
quires  a  vector  rather  than  a  scalar  solution.  Then  in  Sections  III  and  IV  the 
results  are  applied  to  a  turbulent  boundary  layer  on  a  ground  plane  and  on  a 
cylinder. 


*  Note  that  if  the  EM  scattering  process  can  be  fully  described  by  a  single  ray 
in  the  high  frequency  limit,  then  vector  and  scalar  solutions  are  equivalent. 


1 


II.  SCATTERING  FROM  TURBULENT  MEDIA 
IN  THE  GEOMETRICAL  OPTICS  LIMIT 


Limits  on  the  validity  of  vector  solutions  and  scalar  solutions  may  be 
established  in  the  WKB  limit  using  Eq.  (1)  which  includes  both  vector  and 
scalar  solutions  as  special  cases.  Equation  (1)  is  derived  in  the  Appendix. 
For  our  purposes  it  is  sufficient  to  regard  the  constitutive  parameter  K(r)  as 
a  scalar -K(  r  )  =  e(  r  )/e  .  Equation  (1)  is  valid  for  scattering  from  a  weakly 
turbulent  medium  (A  •  In  K(  r )  «  1  and  d/dt  K(  r )  «  d/dt  E(  r  )  in  which  the 
electric  field  E(  r)  is  completely  continuous  except  at  perfectly  conducting 
surfaces. 


<o-<ks  |  lc1)> 


d3"r' 


<6K(  r  )  6K(r' )  E(  r  ) 


•  U(r)  E(r')  •  U(r')>  ,  (1) 

where  <5K(  r  )  =  K(  r )  —  1,  and  where  E  (  r )  is  the  solution  to  the  propagation 

problem  for  a  plane  wave  incident  in  the  direction  on  the  configuration  of 

perfect  conductors  and  turbulent  plasma  media.  U  (  r)  is  the  solution  to  a 

special  propagation  problem  of  a  plane  wave  incident  in  the  — k  direction  on 

s 

the  same  configuration  of  perfect  conductors  but  with  the  plasma  replaced  with 
free-space.  The  function  U  (  r  )  is  essentially  the  Green's  function  of  the  scat¬ 
tering  problem  evaluated  at  r'  —  °°. 

In  the  geometrical  optics  limit,  the  solutions  to  the  propagation  problems 
may  be  expressed  as  a  summation  of  rays. 


(2) 


(3) 
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where  pm  and  pn  are  the  polarization  vectors  of  each  ray  and  the  path  integrals 
are  unique  for  each  ray.  Upon  substitution  of  Eqs.  (2)  and  (3)  into  Eq.  (1),  one 
obtains 


<cr(  k  |  k.  )>  =  2~-  Z  (P  '  P  )  (P  i  *  P  i ) 
N  s'  4tt  u  m  *n  *m'  *n' 


m,m',n,n' 


X 


Em(r)Em'(r)Un(r)Un'(r) 


x  exp 


i(k  -  k 
m  m 


'"V  V'(r+J 


6K( t ) dT 


) 


I 


X  \  d  A  r  6K(  r )  6K(  r  +  A  r  )  exp  { i(  km  —  kn) 


_  n  r+A  r  _  _  p  r 

A  r  +  \  6K(  t  )  dr  —  \  <5K(r)dT 


(4) 


where  A  r  =  r  —  r' ,  and  Em(  r  +  A  r  )  =  Em(  r  )  and  U  (  r  +  A  r  )  =  Un(  r  )  is 
assumed.  In  order  to  relate  Eq.  (4)  to  the  conventional  spectrum  function: 


S(  k  )  =  y  d^  Ar  <<5K(  r  )  6K(  r  + Ar)>  exp{ik  ■  Ar] 


(5) 


it  is  necessary  that 


exp  li(k  —  k  ) 
^  1  '  m  n' 


r+A  r  _  p  r  _ 

SK(t)  dT-  \  6K(r  )  dT 


) 


«  1 


or  conservatively, 


<<5K(7r>  kA  «  2 


(6) 
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where  A  is  the  turbulence  scale  size  and  |  km  —  kn|  ^  k.  Equation  (6)  is 
consistent  with  the  treatment  of  turbulence  as  weak  fluctuations  in  a  laminar 
medium  in  the  Appendix. 

The  two  limiting  cases  of  a  vector  solution  and  a  scalar  solution  may  be 
established  by  further  examination  of  Eq.  (4).  A  vector  solution  is  obtained  if 
the  lengths  of  all  ray  paths  in  the  medium  are  so  short  that  phase  is  not  ran¬ 
domized  by  the  weak  fluctuations.  This  limit  is  given  by 


exp 


+ 


5K( t )  d  r 


«  1 


or  conservatively, 

<<5K(7)2>  kf  «  1  ,  (7) 

where  £  is  the  longest  ray  path  in  the  plasma  medium.  It  is  interesting  to 

compare  this  to  the  conditions  for  validity  of  the  Born  Approximation. 

2 

Salpeter  and  Treiman  indicate  that  the  restriction  on  use  of  the  Born  Approx¬ 
imation  is  usually  the  requirement  that  the  incident  beam  is  not  diffused  over 

z  z  z 

a  wide  angle  <(<50  )>  by  small  angle  scatterings.  Since  <(S0  >  =  <(SK(r)  >  i/A, 
the  Born  Approximation  is  valid  in  the  region 

<6K(~r)2>  £  «  A  .  (8) 

Because  scale  size  A  is  always  small  compared  to  £ ,  Eq.  (8)  is  more  re¬ 
strictive  than  Eq.  (7)  and  the  Born  Approximation  always  requires  a  vector 
solution,  but  the  converse  is  not  necessarily  true.  In  the  limit  given  by  Eq.  (7), 
Eq.  (1)  reduces  to  Yakimenko's  relationship, 


A  scalar  solution  is  obtained  from  Eq.  (1)  in  the  opposite  limit  that  ray 
paths  are  long  and  phase  is  randomized.  In  this  case, 


exp 


i(k  —  k  ,  —  k  +  k  , )  •  \ 
m  m'  n  n'  J 


r  _ 

6K(  r  )  d  T 


»  1 


(10) 


for  all  values  of  m,  m' ,  n,  and  n'  except  m  =  m'  and  n  =  n'.  This  limit  is 
— ►  ? 

obtained  if  </5K(  r )  >  kf  »  1  for  all  ray  paths  where  Eq.  (1)  reduces  to 


<a(ks|  k.)>  =  ^  <6K(T)2>  E  S<Pm  '  P„>2  S<km 

m  n 


k  ) 
n 


x 


rIm(r)  yr) 


(11) 


where  Im(  r  )  =  Em(  r )  E^(  r  ),  and  In(  r  )  =  Un(  r  )  U*(  r  ).  The  scalar  solution 
is  essentially  the  vector  solution  with  cross  terms  removed.  The  ray  inten¬ 
sities  have  been  left  as  functions  of  position  because  diffusion  of  the  incident 
beam  may  be  important  when  <6K(r)  >  kf  »  1  in  a  weakly  turbulent  medium. 

III.  SCATTERING  FROM  TURBULENT  BOUNDARY  LAYERS 
ON  GROUND  PLANES 

The  thicknesses  of  turbulent  boundary  layers  are  normally  sufficiently 

small  for  a  vector  solution  to  be  valid,  except  perhaps  near  grazing  angles  of 

incidence.  Although  the  condition  for  the  validity  of  the  Born  Approximation, 
2 

<6K(r)  >  f  «  A,  is  even  more  restrictive  than  Eq.  (7),  we  choose  to  use  the 

Born  Approximation  for  mathematical  convenience.  In  the  Born  limit,  the 

integrand  of  Eq.  (9)  is  given  by  Eqs.  (2)  and  (3)  with  5K(  r  )  =  0  and  Un(  r )  = 

E  (r)  =  1. 
m 

Horizontal  Polarization  (electric  vector  normal  to  the  plane  of  incidence) 


E  (  r  )  •  U  (  r  )  = 


ik. 


i  k! 


—  e 


-i  k  *  r 
s 


-ik'  •  r 
s 


(13) 
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Vertical  Polarization 


E  (  r )  •  U  (  r  )  =  cos  ©  cos  0 . 

S  1 


+  sin©  sin©, 
s  1 


ik.  -  r  ik'.-  rl  T  -ik  •  r  -ik'  •  r" 
1  ills  s 

i  -  e  e  -  e 


i  k.-  r  i k!-  r 
e  1  +  e  1 


x 


-i  k  •  r  -i  k'  •  r 
e  s  +  e 


where 


/A 


k^  =  k(x  sin©^  -  y  cos©.) 
kj  =  k(x  sin©^  +  y  cos©.) 


k  =—  k(x  sin©  —  y  cos  9  ) 
s  s  s 


k'  =— k(xsin©  +ycos©  )  , 

s  s  s 


and  ©j  and  0g  are  measured  from  the  normal  to  the  ground  plane  y. 
Upon  substitution  into  Eq.  (10)  one  obtains 

Horizontal  Polarization 


"H  =  (??)  <«K2>  Ad  {c/s'  +C1‘S+} 


Vertical  Polarization 


aV  =  (1^)  <6K:2>Ad  [cos2©icos2©s{C1+S-+C1-S  +  } 
+  sin2©,  sin2  0  g  {C^S-  +  C~S  +  ] 

+  sin ©^  sin©g  cos©^^  cos©g  {C^S-  +  Cj  S  +  }j 


(14) 


(15) 


(16) 
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where 


=  s(; 


2k  cos 


9.  ±  0  > 
1  s' 


±  sin  [2kd(cos  0.  ±  cos  0g)]  sin  [2kd  cos  0.] 

^1  -  2  +  kd(cos0.  ±  cos0  )  kd  cos0. 

IS  X 


sin  [2kd  cos  0  ] 


kd  cos© 

s 


±  sin  [2kd(cos  0^  ±  cos  0g)] 

^"2  "  2  +  kd(cos0.  ±  cos0  ) 

1  s 


sin  [ 2kd  cos  0 .] 
kd  cos  0^ 


sin  [2kd  cos  0g] 
kd  cos  0 

s 


±  sin  [2kd(cos  0^  ±  cos  0  g)] 

^  3  "  2  kd(cos0.  ±cos0  ) 

1  s 


and  d  is  the  thickness  of  the  boundary  layer  and  A  is  the  illuminated  surface 
area. 

The  properties  of  this  solution  become  apparent  when  evaluated  in  the 
limit  of  an  optically  thin  boundary  layer  (kd)  «  1  and  compared  to  the  cross 
section  <jQ  of  a  boundary  layer  that  is  not  backed  by  a  ground  plane. 


CTo  =  (?f)  <<5K2>  AdS^2k  cos  9i  2  0-S)  .  (17) 

The  spectrum  of  the  turbulence  is  isotropic  because  the  scale  length  of  the 
turbulence  is  small  compared  to  thickness  of  the  boundary  layer  (i.  e.,  kA  «  1). 
For  the  example  of  a  Kolmogoroff  spectrum, 


S  |2k  cos 


0.  ±  0 
i  s 


=  15.5  A3  | 


1  +  4k2 A2  cos2 


r-4^)] 


-11/3 


~  15.5  A3 
kA  «  1 


(18) 
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For  (kd)2  «  1  and  using  Eqs.  (17)  and  (18),  Eqs.  (15)  and  (16)  reduce  to 


2  2 
16  sin  sin  0 


s 


+  ^  (kd)2  sin©,  sin© 

j  is 


2  ~  2  ^ 
cos  0.  cos  0 


s 


(19) 


!h 

<7  7 

°  (kd)  «1 


16  n  ,A  2  _  2  _ 

T  (kd)  cos  ©.  cos  © 

3  '  l  s 


(20) 


We  see  that  the  effect  of  the  ground  plane  is  virtually  to  eliminate  horizontally 
polarized  scattering  from  thin  boundary  layers  while  enhancing  vertically 


polarized  scattering  for  angles  of  incidence  and  reflection  greater  than  ir/4. 

2 

If  the  medium  is  so  underdense  that  the  <SK(r)  >  l  «  A  condition  is  satisfied 


for  both  i  =  d/cos  0.  and  f  =  d/cos  0  ,  then  vertical  polarization  is  enhanced 
by  12  dB  at  grazing  angles  of  incidence.  If  circularly  polarized  waves  were 
incident  at  grazing  angles  of  incidence,  equal  principal  and  orthogonal  returns 
would  be  scattered  with  a  9-dB  enhancement. 

IV.  SCATTERING  FROM  TURBULENT  BOUNDARY  LAYERS 

ON  PERFECTLY  CONDUCTING  CYLINDERS 

The  incoherent  scattering  from  a  thin  boundary  layer  wrapped  around  a 
large  perfectly  conducting  cylinder  shall  be  determined  in  the  limit  that  ka  »  1 
and  kd  «  1,  where  a  is  the  radius  of  the  cylinder.  In  this  limit,  the  surface 
of  the  cylinder  may  be  assumed  to  be  locally  plane  and  the  results  of  Section  III 
applied.  For  mathematical  convenience  only  backscattering  is  considered. 

At  every  point  on  the  cylinder  the  scattering  properties  of  the  optically  thin 
boundary  layer  may  be  represented  as  a  spectrum  function  as  follows: 


(21) 
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(22) 


2  4 

Srh  =  *6  cos  sin  P  S(2k) 


.  4 


where 


^VH  =  ^HV  =  ^  COs  °H  cos  s*n  P(cos  °h  cos  °V 

+  sino^.  sinOy)  S(2k) 

cos  =  cos  <p  sin© 

cos  =  sin  (p  sin© 

cos  /?  =  cos  (p  cos  © 


(23) 

(24) 

(25) 

(26) 


and  0  is  the  angle  between  the  incident  ray  and  the  normal  to  the  cylinder 
axis  (the  complement  of  the  aspect  angle)  and  <p  is  the  angular  location  on  the 
cylinder. 

The  incoherent  cross  section  is  then  obtained  by  integrating  <p  over  the 
upper  surface  of  the  cylinder,  the  lower  half-surface  being  shadowed. 


cr 


HH  q.2_  ,  2 _  .  20 

-  =  8  sin  0—6  cos  0  sin  0 


o 


'VV 


—  5  sin^0  cos^0(l  —  cos^0) 


2  2  2 
=  8  sin  0+6  cos  0(1  +  sin  0) 


(27) 


,  .  4  2  .  ,  35  2  ~  . 

+  sin  0  cos  0(1  +  -g-  cos  0)  , 


°VH  aHV  ,  .  6rt  ,  1  .  4  n  4n 

-  =  -  =  2  sin  0  +  -=-  sin  0  cos  0 

%  %  2 


o0=  (4J  <«K2>irdbS(2k)  , 
where  b  =  length  of  the  cylinder  illuminated. 


(28) 

(29) 

(30) 
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In  the  cylinder  case,  Oyy  ~  near  grazing  angles  of  incidence  with 

9  dB  of  enhancement  due  to  the  presence  of  the  perfectly  conducting  cylinder 
surface.  The  VV  return  comes  primarily  from  the  top  of  the  cylinder  while 
the  HH  return  comes  primarily  from  the  sides  of  the  cylinder.  Contrary  to 
the  case  of  a  flat  ground  plane,  the  turbulent  boundary  layer  is  polarizing  at 
broadside  incidence. 

Backscattering  from  a  turbulent  boundary  layer  on  a  conical  surface  for 
axisymmetric  incidence  is  similar  to  the  cylinder  case  at  grazing  angles  of 
incidence  except  for  the  lack  of  a  shadow  boundary.  The  point-by-point  spec¬ 
trum  functions  are  given  by  Eqs.  (21),  (22),  and  (23)  but  with  Eqs.  (24),  (25), 
and  (26)  replaced  by 

cos  G!y  =  cos  <p  cos©  cosy  —  sin 0  sin y  ,  (31) 

cos  n =  sin  (p  cos 0  cosy  -  sin©  sin y  ,  (32) 

cosj3  =  cos  <p  sin©  cosy  +  cos©  siny  ,  (33) 

where  y  is  the  one-half  the  vertex  angle  of  the  conical  surface.  After  in¬ 
tegration  over  the  illuminated  surface  of  the  cone,  the  incoherent  cross  section 
is  obtained 


a 


HH 


a 


VV 


O  4  /  b  .  2 

=  8  cos  y  —  6  cos  y  sin  y 


(34) 


rVH  "  °HV 


=  0 


(35) 


Thus  at  grazing  angles  of  incidence,  0  «  (tt/2)  -y,  the  cylinder  solution  to 
zero-order  in  0  equals  the  conical  solution. 


10 


APPENDIX 


Consider  a  plane  wave  E  x(  r )  incident  on  a  weakly  turbulent  medium 
containing  plasma.  The  medium  may  be  bounded  by  and/or  contain  perfectly 
conducting  surfaces.  The  fields  in  the  turbulent  medium  must  satisfy  a  wave 

equation  having  a  variable  or  fluctuating  constitutive  parameter  K(r,  t).  Fol- 

3  — s  ~ 

lowing  Kodis,  an  integral  equation  may  be  set  up  for  the  fields  E  (  r )  scat¬ 
tered  from  the  medium  assuming  weak  fluctuations: 

V  •  ln|Tc(7,t)|  «  1  ,  |Tt(7,t)|  «  ^  [E(7,t)|  , 

esu)  =  J d37T(7 |7)  •  7(7') 


+  (integrals  over  surfaces  of  discontinuous  E(  r))  (A-l) 

where  VXVX  T(r|  r' )  —  k2  T(r|  r')  =  I  <5(r  —  r')  with  n  X  T  =  0  at  the 
conducting  surfaces  defines  the  dyadic  Green's  function,  where 

7=  k2(lt(7)  -  T)  •  E(7)  ,  (A-2) 

and  where  E  (  r )  is  the  total  field  in  the  turbulent  medium. 

The  treatment  of  the  surface  integral  in  (A-l)  is  not  understood.  There¬ 
fore,  the  following  derivation  is  restricted  to  problems  where  sharp  interfaces 
between  the  plasma  and  vacuum  do  not  occur.  The  dyadic  Green's  function 
may  be  determined  by  writing  Eq.  A-l  in  a  form  that  is  equivalent  to  the 
reciprocity  theorem. 

Ps  •  j1  d3  r{Es(r')  6(r-  r')  -  T  (  r  |  r' )  *  J  (  r' )}  .  (A-3) 

Since  E  s(  r )  may  be  considered  to  be  the  field  produced  by  the  source  J  (  r  ) 
in  the  original  configuration  of  plasma  media  and  perfectly  conducting  surfaces 
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but  with  the  plasma  medium  replaced  by  free  space,  then 


Eq(  r  |  7 )  =  pg  •  T  (7  |  7' ) 


(A-4) 


must  be  the  field  produced  by  a  source  p  <5(  r'  —  r )  in  the  same  configuration 
of  conducting  surfaces  and  free-space.  Here  pg  is  the  polarization  vector  of 
the  scattered  fields.  Barring  discontinuous  fields,  Eq.  (A-l)  may  then  be 
written  as 


A 


P 


s 


•  ES(  r) 


J(r') 


(A-5) 


and  the  incoherent  scattering  cross  section  of  the  turbulent  medium  may  im¬ 
mediately  be  obtained  from  Eq.  (A-5) 


<cr(  ks  I  kj)>  =  4tt  |  r  | 2 


ps  •  es(7) 
Pt  •  e\7) 


<c(  kg  |  k.)> 


47T 1  r  |  ^  r 

I  Pi*  E1(7)|'2j 


I d3  r'' 


<Eo(7|7)  •  7(7)  Eo(r  I  7")  •  7 (  r" )> 


(A-6) 


It  is  mathematically  convenient  to  approximate  Eq(  r  |  r’ )  in  the  neighborhood 
of  the  plasma  medium  by  the  solution  U  (  r )  to  the  propagation  problem  for  a 
plane  wave  incident  in  the  direction  —  kg  (from  the  direction  of  the  source  at 
r  as  r  —  00 )  upon  the  free-space  configuration. 


E  (  r  I  r') 
o 


r>j 


r 


— ►<*) 


U(r') 
47r  |  r  | 


(A-7) 
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Then  using  Eq.  A-2  and  assuming  all  incident  fields  have  unit  magnitude,  a 
generalization  of  Yakimenko's  relationship  is  obtained. 


<a(kg|  k.)>  =  J  d3 "r"  j*  d3T'  <[6  K(  r")  •  E(7")  •  U(7")l 


X  [6  K*(  r')  •  E*(r')  •  U*(r')l> 


(A-8) 


where  <5  K(  r' )  =  K(  r' )  -  I  and  E(  r' )  is  the  solution  to  the  propagation  problem 

in  the  configuration  of  plasma  medium  plus  perfectly  conducting  surfaces  for 

incidence  in  the  k^  direction,  and  U  (  r' )  is  the  solution  to  the  propagation 

problem  in  the  same  configuration  but  with  incidence  in  the  —  k  direction  and 

s 

the  plasma  medium  replaced  by  free-space. 
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